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ABSTRACT 

Context. Double-barred galaxies account for almost one third of all barred galaxies, suggesting that secondary stellar bars, which are 
embedded in large-scale primary bars, are long-lived structures. However, up to now it has been hard to self-consistently simulate a 
disc galaxy that sustains two nested stellar bars for longer than a few rotation periods. 

Aims. The dynamical and physical requirements for long-lived triaxiality in the central region of galaxies still need to be clarified. 
Methods. N-body/hydrodynamical simulations including star formation recipes have been performed. Their properties (bar lengths, 
pattern speeds, age of stellar population, and gas content) have been compared with the most recent observational data in order to 
prove that they are representative of double-barred galaxies, even SBO. Overlaps in dynamical resonances and bar modes have been 
looked for using Fourier spectrograms. 

Results. Double-barred galaxies have been successfully simulated with lifetimes as long as 7 Gyr. The stellar and gaseous distri¬ 
butions in the central regions are time dependent and display many observed morphological features (circumnuclear rings, pseudo¬ 
bulges, triaxial bulges, ovals, etc.) typical of barred galaxies, even early-type. The stellar population of the secondary bar is younger 
on average than for the primary large-scale bar. An important feature of these simulations is the absence of any resonance overlap 
for several Gyr. In particular, there is no overlap between the primary bar inner Lindblad resonance and the secondary bar corotation. 
Therefore, mode coupling cannot sustain the secondary bar mode. Star formation is identified here as possibly being responsible for 
bringing energy to the nuclear mode. Star formation is also responsible for limiting the amount of gas in the central region which 
prevents the orbits sustaining the secondary bar from being destroyed. Therefore, the secondary bar can dissolve but reappear after 
s; 1 Gyr as the associated wave is persistent as long as central star formation is active. When star formation is switched off the dynam¬ 
ical perturbation associated with the secondary bar needs several Gyr to fully vanish, although the central morphological signature is 
almost undetectable after 2 Gyr. 

Conclusions. Double-bars can be long-lived in numerical simulations with a gaseous component, even in the absence of overlap of 
resonances or mode coupling, provided that star formation remains active, even moderately, in the central region where the nuclear 
bar lies. 
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1. Introduction 


The early idea that any central triaxial massive component em¬ 
bedded in a stellar disc could be efficient enough to bring gas 
close to the nucleus (e.g. |Kormend^| 1982[ ) has for a long time 
turned the attention of both observers and theoreticians to the 
fruitful framework of “bars within bars” ( |Shlosman et al.|1989|l 
preceded by a few works dealing with triaxial bulges (e.g. ]Louis| 
& Gerhard1|1988^ . At that time, the bars within bars scenario 
was sustained by a few cracial observations (e.g. |de Vaucouleurs] 
1974|[Buta|1986| l. Since then, the observation of small samples 
of double-barred galaxies (e.g. |Buta & Crocker||T993[ |Wozniak| 


.et al.|1995[ l laid the groundwork for larger samples, less biased, 
that has refined the global picture by removing from the list of 
prototypes a few false detections made in the optical and adding 
a lot of new examples mainly detected in the near infrared ( |Er-| 
|win|201 1] for a comprehensive review). 


However, the formation of nuclear bars has remained puz¬ 
zling in spite of the efforts made by several groups to perform 
realistic modelling of these complex systems. Our difficulty in 
creating a standard scenario for the evolution of the central kpc 
is certainly not due to the lack of models that show how the stel- 
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lar material can be assembled to form the inner bar, but rather is 
due to the coexistence of several credible physical mechanisms. 
Most of these mechanisms have been studied mainly through 
numerical simulations (N-body and/or hydrodynamical simu- 


lations) although orbit analysis ( 

Maciejewski & Sparke||1997| 

|20001|Maciejewski & Small|2010 

1 has helped a lot to understand 


the foundation of double-bar dynamics. 


The double-bar formation scenarios can first be divided into 
two major classes: either the large-scale bar forms first, i.e. be¬ 
fore the inner one, or the contrary. The first case is sustained by 
early theoretical works ( [Shlosman et al.|1989|l, simulati ons with 
a gaseous component present (Friedli & Martinet|1993[ Combes 


|1994[|Friedli et al.|1996|l, but also pure collisionless simulations 
( |Rautiainen & Salo|1999[|Curir et al.|2006[ ) making the need for 
a gaseous component questionable if the initial stellar disc is dy¬ 
namically hot and embedded in a massive halo. The second case 
was considered hypothetical until |Debattista & Shen| ( |2007] l cre¬ 
ated such examples by fine tuning the initial conditions of their 
collisionless simulations. One can argue that a third case should 


exist, that is the two bars growing simultaneously ( |Saha & Ma- 
ciejewski||2013 


It seems in all cases that the outer bar forms 


slightly before the inner bar. 
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Whatever the formation scenario is, the main concern with 
inner bars formed in N-body/hydrodynamical simulations is 
that they are short-lived because the central gas concentration 
naturally tends to destroy the bar. However, the lifetime of any 
nuclear bar must be long enough to be compatible with the high 
frequency of double-bars: 30 % of barred galaxies or 20 % of 
all galaxies (|Erwin|[201 l|l. Past simulations that were able to 


produce long-lived double-barred models (e.g. Rautiainen et al. 
2002 [ Debattista & Shen|2007|l are either 2D or purely collision¬ 


less. 

We show here that 3D N-body/hydrodynamical simulations, 
implementing classical star formation recipes, are able to create 
long-lived double-barred galaxies. We identify the dynamical 
and physical processes whereby a double-barred system can sur¬ 
vive on several Gyr. 

This paper is focused on the long-term evolution of double- 
barred galaxies. In Sect. we describe our simulations and in¬ 
troduce postprocessing techniques. The evolution of the cen¬ 
tral region is addressed in detail in Sect.|^ As forthcoming pa¬ 
pers dealing with other aspects of double-barred galaxies will 
be based on this topic, we give a detailed description there. We 
compare the properties of simulations with observational con¬ 
straints in Sect. 1^ Finally, we discuss our results in Sect. and 
conclude in the last section. 


2. Description of the numericai simuiations 

We performed several simulations of disc galaxies varying 
mainly the initial disc scale lengths and the gas contents. For 
each simulation, the conditions for the initial stellar and gaseous 
populations were computed as follows. The initial stellar pop¬ 
ulation was set up to reproduce a typical disc galaxy. The po¬ 
sitions and velocities for Ns particles were drawn from a super¬ 
position of two axisymmetrical [Miyamoto & Nagai| ( |1975| l discs 
of mass respectively Mi and M 2 (cf. Table [IJ). TTie shape of 
Miyamoto & Nagai density distribution depends on the choice 
of two parameters, traditionally called a and b. The [Plummer] 
( |191 l| l sphere {a - 0) and the [Kuzminj ( |1956| ) flat disc (b - 0) 
are the two extreme possible distributions. The use of Miyamoto 
& Nagai distributions enabled us to create initial conditions very 
close to typical disc galaxies. The real radial scale lengths are 
respectively li = ai + b and I 2 - a 2 + b. We chose a scale 
height of b - 0.5 kpc common to all simulations. Thus, with 
our choice of parameters (cf. Table [TJ, the first component can 
be viewed as representing a bulge, the second one a disc, with 
the main advantage that there is no discontinuity in either the 
mass density distribution or the gravitational forces. The initial 
velocity dispersions were computed by solving numerically the 
Jeans equations. The gaseous component is represented by Ng 
particles for a total mass of Mg distributed in a Ig scalelength 
Miyamoto-Nagai disc. 

Some differences in the initial conditions should be noted. 
The disc scalelength of SimL {I 2 - 6.5 kpc, hence the "F" for 
"long") is almost twice that of SimS (3.5 kpc, "S" for "short") for 
the same scaleheight (b = 0.5 kpc). The dynamical timescales of 
the stellar bar development should thus be longer for SimL than 
for SimS. 

The dynamical evolution was computed with a particle-mesh 
N-body code that includes stars, gas, and recipes to simulate star 
formation. The broad outlines of the code are the following. The 
gravitational forces are computed with a particle-mesh method 
using a 3D polar grid with {Nr,N^,Nz) - (60,64,312) active 
cells, leading to a vertical sampling of 50 pc. The smallest ra¬ 
dial cell in the central region is 36 pc large. The hydrodynamics 
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equations are solved using the SPH technique, following closely 
the implementation suggested by |Benz| ( |1990| l. Since we used a 
log-polar grid, we improved the pre-computation of self-forces 
by subdividing each cell in (n^, n^, n^) - (32,6,6) subcells. Self¬ 
forces were then linearly interpolated before being subtracted 
from gravitational forces. The spatial and forces resolutions are 
thus much higher than in our previous studies based on the same 
code (e.g. [Michel-Dansac & Wozniak|2006| Wozniak & Cham- 
pavert|2006l [ 


The star formation and feedback modelling is based on the 
instantaneous star formation approximation (see Michel-Dansac 
|& Wozniak|2004| and reference therein, for details). The major 
steps are 1) the identification of the regions where star formation 
can be ignited, 2) the conversion of a fraction of gas into stars, 
and 3) the computation of the amount of energy and metals in¬ 
jected in the interstellar medium (i.e. the energy and chemical 
feedback from type II supernovae - SNII). Because the last step 
leads to gas heating, a simple treatment of radiative cooling is 
implemented in the energy equation. 

The first task is to identify the gaseous particles that will 
form stars. Friedli & Benz ( 1993| l examined some possible cri¬ 
teria. Not surprisingly, they found that the standard Jeans insta¬ 
bility criterion can be applied to spherical non-rotating gaseous 
systems, but for rotating flat discs that Toomre’s instability crite¬ 
rion ( |Toomre|1964] l is a better indicator. Observational evidence 
also appears to support the use of this criterion as a good in¬ 
dicator for locating star formation at intermediate scalelengths 
( |Kennicutt]|1998| and references therein). In all cases, a parti¬ 
cle i will be assumed to undergo a star formation episode if the 
following condition is verified. 


Qi 


SjKj 


< 4 , 


where Q^. is Toomre’s parameter, i, is the local sound speed, /c, 
is the generalized epicyclic frequency, and F? is the gas surface 
density. The constant A equals unity in the case of an axisymmet- 
ric gaseous disc subject to radial instabilities (cf. Toomre|1964| l. 
However, a value of T 1.4 (with s - 6 km s ) has been de- 
rived from observations ( |Kennicutt|| 1 990) . Since we intend to 
reproduce realistic conditions for star formation, we adopted the 
same value for A. 

The pattern speed Qp of the large-scale bar is easily deter¬ 
mined using its position-angle measured every 10 Myr during 
the simulation run. Any post-processing technique applied on 
snapshots (spaced by 50 Myr) or on movies (timestep of 1 Myr) 
gives the same results. For the nuclear or inner bar, the determi¬ 
nation of Ds is more difficult. The initial size of the nuclear bar 
is unpredictable, and so it is difficult to catch its position-angle 
in real time during the simulation run. The snapshot sampling 
(every 50 Myr) is too long to make a proper measurement as the 
rotation period of the nuclear bar can be much shorter. There¬ 
fore, the nuclear position-angle was determined directly from 
the movies using the inertia moment technique confirmed by 
eye measurements. However, the errors can be large (typically 
2 Myr on rotation period measurements), especially when the 
nuclear bar rotation period is of the order of a few tens of Myr. 

A Fourier transform technique (as used by [Quillen et aLpOlT) 
is unapplicable here because of the frequency cutoff. Therefore, 
for a few periods of interest, simulations have been partly re¬ 
computed using a smaller output timestep (1 Myr) in order to 
apply this technique (see Sect. |^. For these particular periods, 
this technique confirms our manual measurements well within 
our conservative error bars. 
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Table 1. Main initial parameters: name of the run (Model), simulation length, numbers of stellar (Ns) and gas (N^) particles, masses and scale- 
lengths of the two Miyamoto-Nagai initial distributions (Mi,M 2 ,li,h), mass and scalelength of the gas distribution (Afg,/g). The last column 
contains the number of stellar particles at the end of the simulation (Nf"^). 


Model 

End 

Ns 

Ns 

Ml 

M2 

h 

h 

Mg 


^end 


(Myr) 

xlO® 

xlO^ 

xl0“ Mo 

xlO" Mo 

(kpc) 

(kpc) 

xlO" Mo 

(kpc) 

XlO® 

SimL 

9486 

2.5 

5 

0.1 

1.0 

1.5 

6.5 

0.11 

6.0 

3.32 

SimS 

5798 

2.5 

5 

0.1 

1.0 

1.0 

3.5 

0.11 

3.0 

3.21 


The terminology used in the case of double-barred galaxies 
is wide and thus confusing. We thus clarify hereby what terms 
will be used throughout this paper. The simulations will repre¬ 
sent double-barred galaxies, which means two stellar bars will 
coexist. A primary bar is the main bar with respect to its size. It 
is a large-scale structure. It is often named large-scale bar. This 
does not imply that it appears first, even if this is the case for our 
simulations. A secondary bar is a smaller stellar structure em¬ 
bedded in the primary bar. However, it is often called the inner 
or nuclear bar. Hereafter, in the context of this paper, we use the 
word “nuclear” for a rapidly rotating secondary bar extending 
less than 500 pc across, and reserve the term “inner” for larger 
secondary bars. 


3. Evolution of the central regions 

In Figs. and 1^ we show the face-on projections of the central 
region (2x2 kpc) for the two runs used throughout this paper. 
Figure [^focuses on the 1x1 kpc of SimS during the formation 
of the nuclear bar. The stellar mass distribution has been pho¬ 
tometrically calibrated using the technique of 
I Wozniak] ( |2004| l . A blue photometric band (B 
chosen to emphasize the regions with the youngest stellar popu¬ 
lation. For this purpose, the age of the initial stellar population 
has been kept fixed at 10 Gyr for all snapshots. This implies that 
a direct photometric comparison of snapshots has no sense since 
the initial stellar population is not getting old and thus its lumi¬ 
nosity is not dimming. Because it is a post-processing technique, 
this obviously has no impact on the self-consistent chemody- 
namical evolution of the simulation. 

The two simulations show different histories of star forma¬ 
tion (Fig.[2l. Their dynamical history is obviously different, the 
main driver being the dynamical timescale that is roughly dou¬ 
bled for SimL with respect to SimS. 

3.1. Run SimS 

3.1.1. The nuclear bar formation 

The initial disc quickly develops a typical strong bar and a spi¬ 
ral structure both in the stellar and the gaseous components. The 
gravity torques due to the bar and spiral structure drive the gas in¬ 
wards and reorganize the mass distribution even for the old stel¬ 
lar population; this gas inflow occurs in a rather short timescale 
since the star formation rate peaks at f = 0.45 Gyr (cf. Fig.[T]). As 
expected, star formation also occurs along gaseous spiral arms 
that are not shown in Figs. an d0 

>From f ai 0.36 Gyr, the gas distribution in the central 400 pc 
(Fig. 0 starts to twist by an increasing angle that amounts to 
=» 45° at f = 0.4 Gyr. This twist is also visible in the distri¬ 
bution of the new stellar population. At the beginning of the 
twisting process, the new stellar population is aligned with the 
gas distribution. Afterwards, the gas twist angle still increases. 


Michel-Dansac & 
Cousins) has been 



0 2 4 6 8 

Time (Gyr) 


Fig. 1. Star formation rate (SFR) averaged in a sliding window of 
10 Myr for runs SimS (red) and SimL (black). 


quickly reaching 90° at f 0.42 Gyr, but not the stellar angle 
since the collisionless component is less reactive to the torque 
induced by the local gravitational potential twist. Moreover, the 
stellar mass trapped in that region also increases as a result of 
the on-going star formation, thus making the decoupling of the 
stellar component easier. When the gas twist angle continues to 
increase so that the central gas distribution eventually realigns 
with the larger scale flows, the new stellar population is rather 
axisymmetric. 

At f ai 0.5 Gyr there is, however, as much mass in the new 
stellar component as in the gas component. A second episode 
of twisting then begins. At f 0.51 Gyr, the angular mo¬ 
tion of the new stellar population distribution becomes elon¬ 
gated and definitively decouples from the gas distribution. For 
t ^ 0.52 Gyr, the central gas twist is again almost perpendicular 
to the major-axis large-scale stellar bar. At that point, one can 
clearly identify a so-called nuclear secondary bar whose rotating 
motion seems independent from the gas motion. 

During this formation phase, the pattern speed of the nuclear 
bar is very high, leading to a low rotation period of Ts =» 10 Myr 
(Fig.|4|). This value is the dynamical signature of the gas from 
which the new stellar population has been created. The radius of 
the nuclear bar is then /s 250 pc, whereas Ip 4.55 kpc, which 
leads to l^/lp ~ 0.056. 
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Fig. 2. Run SimS: gas mass surface density (left) and B-Cousins surface brightness (right) in the central ±1 kpc for t < 0.7 Gyr. Isocontours 
are spaced by 0.15 magnitude in logMopc^^ for the gas and 0.5 Bmagpc^^ for the stars. The map resolution is 10 pc. Gas particles have been 
convolved by the SPH kernel to obtain the real distribution. Gas surface densities below 1 Mq pc“^ are contoured with dashed lines. The straight 
dotted line is the primary bar position-angle. 


3.1.2. The nuclear bar evolution 

Once the nuclear bar has been formed and rotates independently 
of the gas (what was called “decoupling” in some studies from 
the 1990s), the central gas concentration quickly gets aligned 
with the stellar component. For t ^ 0.63 Gyr the gas dynamics 
can be considered fully driven by the nuclear bar evolution (ori¬ 
entation, velocity) since the gas no longer represents the most 
massive component in that part of the galaxy. The morphol¬ 
ogy is modified consequently. Indeed, the gas distribution is 
more ring-like around the nuclear bar or sometimes disc-like. 
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especially when the energy feedback from SNII temporarily dis¬ 
solves a part of the gaseous ring. 

After 1 Gyr, the total mass inside the central kpc has in¬ 
creased by a factor of 1.5. Indeed, a nuclear gas disc is formed 
from the accumulation of gas in the centre, and new stars are 
actively formed there. The first cause of this mass inflow is the 
overall reorganization of the mass distribution under the influ¬ 
ence of the stellar bar, even the old population. Owing to grav¬ 
itational torques exerted on the gas by the stellar bar, the extra 
mass in the form of gas and new stars amounts to 1.97 10^ Mq 
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Fig. 3. As in Fig.[^ but for t > 0.7 Gyr in the central ±2 kpc. The map resolution is now 20 pc. Gas isocontours are spaced by 0.25 magnitude 
in log Mq pc“^. The radius of ILRp is indicated by dashed white circles. 


at f = 1 Gyr, which is only 31% of the whole additional mass. 
The redistribution of the old stellar population contributes to the 
other 69%. 

The length of the nuclear bar slowly increases as it slows 
down. The large-scale bar also increases, but at a lower rate. 
Therefore, at f = 2 Gyr k ~ I kpc, /p 5 kpc so that l^/l^ ^ 0.2. 

3.1.3. The nuclear bar dissolution. Formation of the nuclear 
disc 

While the nuclear stellar bar continuously grows, the gaseous 
ring radius increases. For f ^ 2 Gyr the nuclear bar gets also 


thicker, looking sometimes like an oval surrounded by a ring 
of gas. The ring itself gets broader as the nuclear bar thickens. 
The nuclear bar is sometimes slightly off-centred, in particular 
around f 1.6 Gyr, but this is a transient effect. 

It is impossible to precisely date the moment when the nu¬ 
clear bar can be considered as dissolved because the dissolu¬ 
tion process is smooth and has a long timescale. Except for 
t ^ 2.7 - 3.0 Gyr, one can consider that the nuclear bar has 
been replaced by a nuclear disc. This stellar nuclear disc has a 
clear signature in the edge-on mass distribution. It also drives 
the shape of the gaseous nuclear ring. Any slight temporary 
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Fig. 4. Period of rotation for the primary bar (in black) and the nuclear Fig. 6. As in Fig.|^for SimL. 
bar (in red) for SimS. 


ovalization of the stellar nuclear disc leads to an alignment of 
the major-axis of the gaseous ring. 

Since the gas accumulation in the central region continues, 
whereas the star formation rate is very low in the ring, the gas 
mass increases in the ring and the nucleus. For f > 4 Gyr the 
gaseous ring becomes unstable. It starts to twist and collapse. 
This permits the gas flows to reach the nucleus and leads the star 
formation rate to briefly increase (Fig. [TJ. As a result, the nice 
broad gaseous ring is replaced by a gaseous disc of the same size 
as the nuclear disc for t > 4.3 Gyr until the end of the simulation. 

3.2. Run SimL 

3.2.1. The nuclear bar formation: the first secondary bar 

As for SimS, the initial disc of SimL develops a typical large- 
scale strong bar and spiral structure. However, because the disc 
scalelength of SimL is roughly twice that of SimS, the dynam¬ 
ical timescale for the development of the main stellar bar is 
significantly greater. It is only after t ^ 0.8 Gyr that the new 
stellar population created in the central gas concentration forms 
a nuclear bar-like structure (Fig. [^. This nuclear structure is 
very small, ^ 250 pc, whereas /p ^ 5 kpc, which leads to 
k/lp = 0.05. It is also transient, lasting less than =» 0.6 Gyr, 
and is often asymmetric and/or off-centred. It then dissolves 
into a stellar nuclear spiral structure surrounded by a circum- 
nuclear gaseous ring. It must be also mentioned that the large- 
scale bar growth is not completed yet. However, even if the 
length of both bars increases during this phase, the ratio 4/^p 
remains approximately constant. For instance, at t ^ 1.3 Gyr, 
4//p ~ 0.3/6.0 = 0.05. 

Because of its time dependent morphology, the pattern speed 
of this nuclear structure is difficult to measure precisely. The pat¬ 
tern rotation period of the large-scale bar and the nuclear struc¬ 
ture are shown in Fig.[^ It was not possible to determine prop¬ 
erly the nuclear rotation period for the whole lifetime because of 
the quick changes in morphology. However, the nuclear struc- 
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ture has a very short rotation period (between 10 and 30 Myr) 
that strongly contrasts with the primary bar rotation period; it 
is at least ten times longer. During the lifetime of the nuclear 
structure, the large-scale bar temporarily slows down, the rota¬ 
tion period increasing from 200 Myr to 300 Myr (Fig.[^. 

When the nuclear structure dissolves, both a gaseous disc 
and a stellar disc (mainly made of new stellar population), sur¬ 
rounded by a ring-like structure, survive. These discs form the 
remnant of the nuclear bar. The size of both gaseous and stel¬ 
lar discs increases to « 1.8 kpc in diameter, which leads to 
Is/lp ~ 0.9/8 = 0.11. However, because the remnant stellar 
disc is no longer barred, l^/lp has another meaning here. Tran¬ 
sient gaseous and/or stellar spiral arms often appear inside and 
outside the ring-like structure until f ^ 2.6 Gyr. 


3.2.2. The inner bar formation: the second secondary bar 


After t ^ 2.6 Gyr, the inner region returns to a bar-like shape. 
The pattern speed can be determined much more precisely now 
than for the first nuclear bar. The gaseous and stellar compo¬ 
nents both display the same features. The inner bar ends in a 
broad circumnuclear ring. In the centre, the remnant of the first 
nuclear bar settles in the centre as a flattened small bulge. As the 
stellar inner bar begins to turn into an oval, the gaseous coun¬ 
terpart progressively dissolves. Only a gaseous circumnuclear 
ring survives, that gradually gets broader (up to roughly 500 pc 
across the torus). 

For f ^ 4 Gyr, it becomes difficult to unambiguously de¬ 
termine a nuclear bar structure. The central kpc region looks 


like the triaxial pseudobulges described by Kormendy & Ken- 
|nicutt| ( |2()04] l. However, even if the contrast between the inner 
structure and the old stellar population background is too low to 
permit a proper pattern speed determination, the inner region is 
not devoid of non-axisymmetries. Therefore, for f ^ 4.5 Gyr, 
the whole inner stellar structure, including the now vanishing 
circumnuclear ring, turns into a strong oval that looks like a 
Is = 2.2 kpc long bar {k/lp ~ 0.275). The shape of this second 
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(in the time sequence) inner bar depends strongly on its relative 
orientation with respect to the large-scale bar. The higher axis 
ratio is obtained when the angle between the two bars reaches 
90°. When the two bars are aligned, the axis ratio is roughly 
halved. 


For t ^ 4.5 Gyr, the responsive gaseous ring also starts to 
turn into an oval, meaning that the gas ring shrinks and expands 
depending on the relative orientation of the inner bar with respect 
to the large-scale bar. The ring progressively dissolves into a 
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disc that closely follows the same changes in morphology as the 
stellar nuclear bar and thus shares the same orientation. 

Because star formation, which mainly occurs in the central 
region, consumes the gas, the nuclear disc size signihcantly de¬ 
creases for f ^ 6 Gyr and becomes dynamically unstable. For 
6 < f < 7 Gyr, the star formation rate temporarily increases 
(Fig. 0 because of the collapse of the gaseous nuclear ring. This 
is the last episode of star formation in that region. At f « 7 Gyr, 
the gaseous nuclear disc has almost fully disappeared, whereas 
the stellar counterpart is made of an inner bar with 4 ~ 1-3 kpc 
embedded in a slightly larger stellar disc. The large-scale bar has 
Ip = 8.5 kpc leading to 4//p =» 0.15. 

As the stellar inner bar (oval-like) still exists for t > 1 Gyr, 
the gas inflow continues in such a way that a gaseous circumnu- 
clear ring forms again. The gas mass surface density in the ring, 
and in the region encircled by the ring, is higher than 1 Mq pc“^ 
, but never reaches the same high values as during the first 2 Gyr. 
In particular, the gas distribution is smoother. 

These structures, the inner stellar bar or oval and gaseous 
circumnuclear ring, last until the end of the simulation at f ^ 
9.5 Gyr without any significant morphological change. The only 
noticeable fact is the permanent slow down of both bars due to 
the classical angular momentum exchange with the large-scale 
disc. 

At the end of the simulation, 4/4 ~ 1-5/9 = 0.17, quite 
close to the value alt - 1 Gyr. Indeed, the two bars slow down, 
but at different rates. 



Radius (kpc) 

Fig. 7. Radial profile of the age difference with the time of the snap¬ 
shot (t = 9.486 Gyr) for SimL. The mean age has been computed in 
concentric annuli of 100 pc wide. For the sake of comparison, the ra¬ 
dius of the NGC 5850 inner bar and the limit of age measurement jd^ 
|Lorenzo-Caceres et al.|201^ have been scaled to SimL and overplotted 
as dashed and dot-dashed lines, respectively. 


4. Comparison with observations 

4.1. Gas mass in the central kpc 


It has been argued in the past ( |Petitpas & Wilson|2004[ and ref¬ 
erences therein) that large amounts of gas are not required to 
form and sustain a double-bar. However, if the nuclear bar is 
long-lived, and we can imagine that the galaxy has been gas rich 
in the past, then its gas has been consumed and possibly redis¬ 
tributed thus explaining the current low content. Alternatively, 
if the nuclear bar is a short-lived phenomenon, one has to define 
the lower gas mass limit to initiate the formation process. This 
means that the expressions ‘gas rich’ or ‘gas poor’ used in the 
literature must be clarified. 

The last development stage (i.e. after 7 Gyr) of SimL is 
morphologically close to SBO galaxies, and so it is morphologi¬ 
cally comparable to the |Petitpas & Wilson| ( |2004) l observations. 
We can thus determine the amount of gas in a region of sim¬ 
ilar surface limited by the beam size. For ten objects, the ra¬ 
dius encircling the mass that has been measured ranges from 
0.3 kpc to 2.6 kpc, but for their SBO subsample (NGC 2859, 
NGC 2950, NGC 3081, NGC 4340, and NGC 4371) r < 1 kpc 
and r < 2.3 kpc are more representative. For r < 1 kpc, observa¬ 
tional gas mass upper limits range from 1.5 10^ to 8.1 10^ Mq , 
whereas it amounts to 4.4 10^ Mq for SimL. For r < 2.3 kpc it in¬ 
creases to 2.8 10^ Mq for our simulation, whereas for NGC 3081 
the gas mass amounts to 6.2 or 6.5 10* Mq. Since our simulation 
is not fine tuned to fit any particular galaxy, we can conclude 
that SimL is well within the observational con straints. A value 
of 4.4 10^ Mq is low enough to explain why Petitpas & Wil 


|son[ ( |2004| l were unable to detect molecular gas in NGC 2859, 
NGC 4340, and NGC 4371, but obviously one cannot conclude 
that the gas component do not play a role in the formation of 
nuclear bars. 


4.2. Stellar population in nuclear bars 

Studying the stellar populations of double-barred galaxies 
NGC 2859, NGC 394 1, NGC 4725, and NGC 5850, |de Loren^ 
[Caceres et al.| ( |2013) have recently concluded that the inner bars 
are younger and more metal rich than the large-scale primary 
bars. It has been argued ( |Wozniak|2007[ ) that the stellar age dis¬ 
tribution in bars has to be interpreted with a dynamical approach. 
For instance, two regions of apparent low age at the end of the 
large-scale bars ( |Perez et al.|2007[ ) are due to the accumulation of 
a composite stellar population, younger in average and trapped 
in orbits shaped like an ellipse in average, aligned with the bar. 

To make a comparison with |de Lorenzo-Caceres et al7| 
( |2013 I, a mean age profile has been computed as the average of 
individual particle age in annuli 100 pc wide and +100 pc thick 
centred on the nucleus. Being an azimuthal average approach, 
each annulus mixes the population inside and outside the struc¬ 
tures (inner and outer bars, circumnuclear ring). Moreover, all 
particles of the initial population arbitrarily have the same age (0 
at the beginning of the simulation) so that snapshot time is also 
the age of the initial population. 

In Fig. 1^ we show the relative mean age of the particle pop¬ 
ulation at f = 9.486 Gyr. The region inside the nuclear bar is 
2.5-3.0 Gyr younger than the large-scale bar. Because of the 
mixing with older regions outside the inner bar, and the rough 
assumption that all initial particles have an age of 9.486 Gyr, 
this value is certainly a lower limit. The outermost region of the 
large-scale bar is also younger. This effect is obviouly that dis¬ 
covered by |Perez et al.| ( |2007) l and explained by |Wozniak| ( |2007| l. 

Even if a direct and quantitative comparison with 
Lorenzo-Caceres et al. p013|l is not possible because the av¬ 


eraging technique is different (they measured ages along ellipses 
fitted on the isochrone maps), the results are in good qualita- 
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Fig. 8. Ratio of the primary to secondary bar pattern speed. Black 
dots are for SimL, red dots for SimS. 


tive agreement. In particular, the age difference for NGC 5850 
amounts to « 4 Gyr between the nuclear bar and the outermost 
measurements. To help the comparison, the large-scale bar ra¬ 
dius of NGC 5850 (using the de Lorenzo-Caceres et al. ( |20131 l 
value, i.e. 63") has been scaled to SimL (/p ^ 9 kpc). The large- 
scale bar radius was used as its determination is much more ac¬ 
curate than the nucle ar radius. In Fig. [7] the oute rmost limit of 
age measurement by |de Lorenzo-Caceres et al.| ( |2013[ l is also 
plotted showing that their age measurements might be well in¬ 
side the large-scale bar. 


4.3. Pattern speeds 


The ratio of pattern speeds (Gp/Qs or equivalently TJTp, cf. 
Fig[^ is not constant over time. For the two simulations, both 
Qp and Qs decrease as a function of time (Figs.fflandl^ but are 
not locked into a particular ratio. The slopes seem different be¬ 
tween SimS and SimL, but are in fact consistent with a factor of 
roughly two in dynamical timescales. This means that the ratio is 
the same for SimL after twice the time needed by SimS. Because 
on average decreases much faster than Qp, the ratio Qp/Qs 
monotonically increases up to 0.21 - 0.25. Only a change in 
the slope and a lower dispersion between 6 and 8 Gyr are notable 
for SimL, mainly due to Qj variations during the circumnuclear 
ring collapse (cf. Sect |3.2[ ). The fluctuations of the ratio is due to 
oscillations of Qj and Qp (Figs.fflandl^ that depend on the rela¬ 
tive phase of the two bars, a fact alrea^ observed by |Englmaier| 
|& Shlosman| ( |2004| and |Debattista & Shen| ( |2007| l. 

The best model of|Friedli & Martinet|(|1993|l has Qp/Qj ^ 
0.32, and 0.3 for |Fnglmaier & Shlosman| ( |2004p , values that our 
simulations never reach. Our values of Qp and Qs are compara¬ 
ble, however, to the few measurements in double-barred galax¬ 
ies. Most observational methods to directly estimate Qp cannot 
be easily applied to Qj and, more generally, are not suited when 
multiple patterns are present (|Maciejewski|200^|Shen & Debat-| 


|tista|200^ . However, keeping in mind the questionable viability 
of Qi measurement, we must rely on available attempts. 

Therefore, if we assume two b ars in NGC 5248 (cf. dis cus- 
1201311, Qn = 30 kms 'kpc 


Sion m 


van der Laan et al. 


Go¬ 


gee et al. 20021 and Qs^ 135 kms 'kpc ' leading to a ratio 

Qp/. 

0.13 

12s ^ 0.23. 

7!orsini et al. (2003 1 found a value as low as 

tor using the Iremame-Weinberg method (but 

see 

Maciejewski 

2006 Shen & Debattista 2009jl. Hernandez 

et a 

1. (200511 

"athi et aLJ2009]l applied the same technique 


to Ha velocity fields to determine pattern speeds and found 
0.3 < Qp/Qj < 0.55, but with large error bars on Qj. An in- 
depth analysis of NGC 1068, based both on multiwavelength ob¬ 
servations and numerical simulations, led |Emsellem et al!] ( |20061 l 
to conclude that 0.19 ^ ^n/^s ^ 0.26. 


Zhang & Buta (20071 introduced a potential-density phase- 


shift method that is able to determine the position of multi¬ 
ple corotations in a galaxy without any kinematical informa¬ 
tion. This method is thus particularly well suited to double- 
barred galaxies. Moreover, it is based on the properties of quasi¬ 
stationary density wave modes. [Zhang & Buta| ( j2007| l did not 
publish values for Qp and Q^, but only radii for CRs and CRp. 
However, the ratio of corotation radii CRs/CRp can be used as it 
evolves like Qp/Qj. Eor their small sample (seven objects), they 
found ratios between 0.086 (NGC 1530) and 0.3 (NGC 936) and 
a mean of =» 0.17. Eor SimS and SimL, CRs/CRp starts at a; 0.03 
when a nuclear bar is just detected for the first time and then 
increases to respectively 0.24 and 0.30. Our values are thus 
comparable even for SBO galaxies. 

The most recent survey ( jEont et al.|2014b| l of pattern speeds 
in double-barred galaxies, using a phase reversal technique in 
Ha kinematical maps ( jEont et al.|2014a| l found a narrow range 
of ratios 0G9/^, a bit higher than our end values. However, 
jEont et aL]p014b| l have introduced an important assumption: the 
nuclear bar ends near its corotation (CRs) which in turn overlaps 
(or is close to) the primary bar ILR (ILRp). If this new tech¬ 
nique could be useful to determine the location of resonances 
(although a firm physical foundation must be found), it gives no 
information on the type of these resonances, especially the main 
ones, ILR, ultraharmonic (UHR, also called 4:1), CR, and OLR. 
Therefore, the identification of the resonance needs to assume 
that any bar ends inside and close to its corotation. Whereas 
this is reasonable true for the primary bars (although it has been 
argued by Michel-Dansac & Wozniak ( 2006| l that typical obser¬ 
vational bar length measurement techniques are more correlated 
with the UHR resonance than with the CR), this assumption is 
questionable for nuclear bars. Indeed, if CRs is determined in 
SimL as being the end of the nuclear bar, and moreover iden¬ 
tical to ILRp, this leads to a large underestimation of Qj. Eor 
the sake of illustration, we consider the moment t - 1.5 Gyr. 
At this time, Rilrp ~ 4.84 kpc and Qp « 15 kms^'kpc"'. If 
we assume that this radius is identical to the corotation radius 
of the nuclear bar, this leads to Q, 45 km s ' kpc ' instead of 


73 km s ' kpc ' (see also Eig. 10 in the next section). There¬ 


fore, the ratio Qp/Qs might be spuriously estimated to be 0.33, 


a value which appears to be in good agreement with Eont et al. 
( j2014bj l, instead of a; 0.2, its real value. 


4.4. Bar lengths 

The absolute sizes of simulated secondary bars fit quite well with 
the values given by Erwin (201 Ijl. Erwin showed that k/lp, the 
secondary over prirnary lengths, ranges from 0.025 to 0.25 
with a median of ^ 0.12. Eor the first nuclear bar of SimL, 
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4//p remains roughly constant at 0.05. It is only when the nu¬ 
clear bar has been fully dissolved that this ratio reaches 0.11 
(this is then the size of the remnant nuclear disc relative to the 
large-scale bar). The inner bar restarts with a rather high ratio 
(0.275 at f ai 4.5 Gyr), but then decreases to 0.15-0.19. How¬ 
ever, the observational relative size distribution does not show 
any strong evidence of two different kinds of secondary bars that 
could help to distinguish the first generation from the subsequent 
ones. Our simulations should also show some similarities with 
that of |Friedli & Martin^ ( |1993| l since our recipes are quite sim¬ 
ilar, although our initial conditions are different. Indeed, the ba¬ 
sic process leading to the formation of a nuclear bar is the same. 
However, their best model has a bar length ratio of 0.26. 

Although it is incorrect to base an argument on absolute 
sizes, it is also noteworthy that, all other factors being equal, 
if primary bars have low Qp (as for SimL), they are longer than 
bars with high Qp (as for SimS), assuming here that they sys¬ 
tematically fill the region inside their corotation radius. [Erwin 
( |201 l| l emphasizes that the sizes of large-scale bars are longer 
in double-barred galaxies than in single-barred ones whereas all 
the other properties are similar. 

5. Discussion 

SimS and SimL having morphological, kinematical and stellar 
population properties similar to real galaxies, we can now dis¬ 
cuss their dynamical structure and long-term evolution. 


5.1. Resonance overlap or not? 


In the framework of the epicyclic approximation, we solved the 
equations Q.- Kfri - Q.p s for the resonance radii (ILRp and ILRj 
for n - 2, UHRp and UHRs f or n = 4, CRp and CRs for n = oo, 
OLRp and OLRs for n - -2). Athanassoula ( 1992|l and Michel- 
[Dansac & Wozniakj ( |2006| l have shown that, even in the case of 
strong bars for which the epicyclic approximation breaks down, 
the errors on the resonance positions remain within 10%, espe¬ 
cially in the case of the ILR and the CR. In Fig. we show the 
evolution of ILR, UHR, CR, and OLR as a function of time. 
Since most of the time two ILRs exist, hereafter we only deal 
with the outermost one (named oILR by most authors). More¬ 
over, because it is very close to the centre in the case of SimS, 
the position of ILRs is not always reliable for f 2 Gyr. 

The immediate observation is that the whole resonance sys¬ 
tem of the inner bar, up to OLRs, remains well inside the ILRp 
for a long time (roughly 2 Gyr for SimS and 8 Gyr for SimL). 
Therefore, there is no CRs-ILRp overlap, unlike Friedli & Mar- 


tinet||(T993| and many other similar studies that followed. 


However, such a situation is not new as a few examples have 
been already reported by [Rautiainen & Salo] ( |1999| l and [Rauti-j 
ainen & Salo ( 2000|l, but with different experimental setups (ei¬ 


ther 2D simulations or an analytical bulge component). More¬ 
over, [Moellenhoff et ah] ( |1995| l reported evidence of an OLRs- 
ILRp coupling in the case of NGC 4736, before the paradigm of 
CRs-ILRp coupling was firmly established. They claimed this 
configuration could weaken the interaction and increase the life¬ 
time of the nuclear bar. 

To further investigate the dynamical impact of possible res¬ 
onance overlaps, we have built spectrograms, as introduced by 
Sellwood & Athanassoula (19861 and widely used (see e.g. 
Quillen et al.|2011| and references therein). We note m the az¬ 


imuthal wavenumber of a wave of frequency oj. Spectrograms 
can be used to determine the pattern speed Q of any wave as be¬ 
ing Q = w/m and helps to identify any potential non-linear mode 
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Fig. 9. Radius of the linear resonances (ILR, UHR, CR, OLR) for 
the secondary bar (subscript ’s’, red triangles) and the primary bar (’p’, 
black dots) in the simulations SimS (top) and SimL (bottom). OLRp is 
not shown. 


coupling when resonances overlap. As discussed by [Rautiain^ 
|& Saloj ( |1999| l and [Quillen et al.j ( [201 1[ ), spectrograms show 
many small features, often transient, probably linked to the main 
modes. Some of these features with a well-defined pattern speed 
can be beat modes triggered by the non-linear couplings of other 


waves at the location of resonance overlaps. Rautiainen & Salo 
\\999\ found many kinds of resonance overlaps in addition to 


p, although they were not always able to 


the classical CRj-ILR, 
determine whether these couplings were stationary or acciden¬ 
tal. 
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Fig. 10. m = 2 power spectra for SiitiL in log scale. The time window 
spans 269 Myr centred at f = 3.3 and t = 7.5 Gyr. The left vertical 
scale gives values of a>, the wave frequency in Myr^*, whereas the right 
vertical scale is for the pattern speed O = tu/2 in km s“' kpc“*. The hor¬ 
izontal white lines represent and Dp, respectively, at the beginning 
of the time window (full line) and the end (dot-dashed one), fl, de¬ 
creases from 146 to 139 at f = 3.3 Gyr and from 81 to 71 km s“' kpc“* 
at f = 7.5 Gyr, whereas Gp decreases by less than 1 kms^' kpc“* so 
that the two lines cannot be distinguished. The averaged curves 212 - k 
(which allows the ILR to be identified) and 2fl + k (for the OLR) are 
drawn as white short dashed lines, 212 - /c/2 (for the UHR) as long 
dashed line, and 212 as a solid line (for the CR). A logarithmic scale has 
been chosen for the radius to emphasize the central kpc. 


Figure[T^shows spectrograms for m - 2, centred at f = 3 and 
7.5 Gyr for SimL and computed in a window of 269 Myr corre¬ 
sponding to 256 outputs. To compute the spectra we included 
only the stellar component to extract the information on stel¬ 
lar modes. Similar spectrograms have been obtained for SimS. 
Spectra with greater wavenumbers (m > 4) confirm what can be 
inferred from m - 2, in particular the presence of harmonics of 
the large-scale bar mode. They are not shown here. For the two 
moments shown in Fig. [^representative of the rest of the sim¬ 
ulations, the two dominant modes correspond to the secondary 
and primary bars, with pattern speed values in good agreement 
with all other measurements (cf. Sect. j^. The features are not 
as thin as for jMasset & Tagg"^ ( \991\ . This is due in part to our 
narrower time window. This also shows that these features do 
not correspond to quasi-stationary structures. 


Figure [^clearly shows that most of the power of the nuclear 
bar mode remains inside the UlTRs. For the large-scale primary 
bar, it has been argued that the most appropriate orbits to sustain 
a bar are those existing inside the UHR (see e.g. Sect. 6.1 of 
[Michel-Dansac & Wozniak|2006| and reference therein). This 
also seems to be true for the secondary bars of SimL and SimS. 
This additionally supports our previous claim (Sect. |4.3] l that one 
cannot systematically assume that the morphological end of a 
nuclear bar coincides with CRj. 

Finally, whatever the approach, the linear resonance analysis 
or the mode analysis both lead to the same conclusion: there is 
no resonance overlap or mode coupling in SimL and SimS. 


5.2. Double-bar lifetime 


|Rautiainen et al.| ( [2002| l reported the case of a long-lived double- 
barred system (their model I) which could be very similar to 
SimL, although SimL is 3D and fully self-consistent whereas 
their model I is 2D and the bulge is represented by an analyt¬ 
ical Plummer sphere. However, the fact that long-lived systems 
can exist is important in order to explain the high frequency of 
double-barred galaxies. In our case, the nuclear bar lifetime is 
clearly linked to the lifetime of the central bar mode. The fact 
that we are not able to determine a proper pattern speed during 
roughly 1 Gyr centred on f « 2 and 4 Gyr (cf Fig. [^ does 
not mean that dynamically the nuclear mode has disappeared. 
Indeed, the mode still exists (Fig. ED at a pattern speed of re¬ 
spectively ^ 215 and 125 kms“^ kpc“^. This means that even 
if, on photometrical grounds, it can be stated that the nuclear 
bar has been dissolved into a triaxial bulge or pseudobulge (j^ 
|rmendy & Kennicutt|2004[ ), or even a weaker structure like nu- 
clear spirals (cf. Sect. [^, the bar mode does survive, which in 
turn facilitates the morphological revival of the nuclear bar later 
on. 


Debattista & Shen (|2007 1 also reported long-lived double- 


barred systems in purely collisionless simulations. A direct com¬ 
parison is not possible since the initial conditions of their sim¬ 
ulations, especially the bulge part, were fine-tuned to produce 
such gas-free double-barred galaxies. However, their enforced 
rotating bulge mimics perfectly the effect of a dissolved or weak 
first nuclear bar. 


Another crucial dynamical feature of our models that could 
play an important role in the lifetime of the secondary bar, is the 
absence of any overlap of resonances for a long time. Indeed, 
the whole nuclear dynamical system evolves as if embedded in 
a quasi-stationary gravitational background due to the rest of the 
galaxy. The mass (both gas and stars) encircled by the OLRs 
is so large that the two dynamical systems are unable to couple 
together through their resonances. For SimL, almost 38 % of the 
total mass of the galaxy lie inside R % 5 kpc -1 Gyr. 

All previous works that have looked for mode coupling in 3D 
self-consistent N-body/hydrodynamical simulations have found 
an CRs-ILRp overlap whereas their double-barred systems were 
short-lived (e.g. Friedli &Martinet|1993| l. One noticeable excep¬ 
tion is |Rautiainen etal. p002| l who argue that this mode coupling 
is not necessary for the coexistence of bars and show examples 
of other kinds of coupling. [Rautiainen & Saloj ( |2000| l, using 
2D N-body simulations and inelastically colliding massless gas 
particles, exhibit an example very close to SimL with OLRs well 
inside ILRp (their model A2.5). This is the situation of SimS 
for t % 2 Gyr and SimL for f :$ 8 Gyr. However, their model 
develops a very weak primary bar (more oval than SimL) and the 
m = 2 of the nuclear mode remains close to the maximum of 
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Fig. 11. m = 2 power spectra for SimL centred at r = 2 (top panel) 
and t = 4.275 Gyr (bottom) in a window of 538 Myr, whereas the inner 
bar cannot be detected on density maps. Consequently, cannot be 
measured because the position-angle of the central structure cannot be 
properly determined. 


Q - a '/2 all through the simulation, which is not the case of our 
simulations. 

Can we infer a strong dynamical perturbation when OLRj 
position comes close to ILRp and eventually overlaps with ILRp? 
For SimS, less than 1 Gyr after OLRj intersects ILRp, the nuclear 
bar is transformed into a disc with a very time-dependent shape. 
It is also remarkable that the nuclear bar then starts to dissolve 
in the case of SimS whereas it continues to evolve as a large oval 
for SimL, showing larger fluctuations in Q.^. Three spectrograms 
centred at f = 2,3.4, and 5.47 Gyr for SimS (Fig.[T^ show that 
the nuclear bar mode progressively vanishes whereas the large- 
scale bar mode (with a lower pattern speed) expands toward the 
centre. However, the OLRj-ILRp overlap seems to play no role. 
Indeed, the primary bar mode is at a maximum just inside ILRj 
, but we have not found a mode associated with the nuclear bar 
that might be present near the OLRj to explain the energy and 
angular momentum exchanges. The overlap thus seems neutral 
and the nuclear bar mode extinction could simply be due to clas¬ 
sical Landau damping at ILRs. 

In the case of SimL, the shape of the inner bar is oval for 
f ^ 7 Gyr. Looking at the spectrograms (Fig. 13 i for t - 7.6 
and 9.2 Gyr, it appears that the power of the inner bar mode is 
concentrated just outside ILRj whereas the maximum power of 
the primary bar is just inside ILRp, as it is for SimS. Again, it 
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Fig. 12. m = 2 power spectra for SimS centred at ? = 2 (top panel), 
when OLRj intersects ILRp, f = 3.4 (middle), and t = 5.47 Gyr (bottom) 
where the nuclear bar is dissolved. For t -2 Gyr the time window has 
been halved (269 Myr) because fl, evolves rapidly. 


is hard to find any mode between CRs and OLRs. Unlike SimS, 
several modes can be detected at low level with < G < or 
Q > Qj, but only inside a region limited by Q(r). As we do for 
SimS, we can suspect Landau damping at ILRs to be responsible 
for the inner bar mode decrease. 
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Fig. 13. m - 2 power spectra for SiitiL centred at f = 7.6 (top 
panel) just before OLRj intersects ILRp, 1 = 9.2 Gyr (bottom) whereas 
the nuclear bar has been replaced by a large rotation oval. The time 
windows is 568 Myr wide. 


5.3. What makes double-bars long-lived? 

The lack of mode coupling is a major difference with previous 
simulations. For pure N-body simulations, but also for models 
with massless gas particles, this situation would prevent the nu¬ 
clear bar from being long-lived. Indeed, apart form the swing 
amplification at corotation that does not operate here, there is no 
mechanism to bring energy to the inner bar wave. Therefore, this 
wave should vanish by Landau damping at ILRj. 

In our simulations, the first difference is the presence of a 
self-gravitating gas component. This component behaves dif¬ 
ferently from the stellar component because of its dissipative 
nature. However, as shown for example by |Friedli & Martinet| 
P993| l, the gas does not help to sustain a long-lived nuclear 
bar because the mass accumulation close to the centre reinforces 
ILRp , which in turn destabilizes the nuclear bar orbits. 

Another major difference is the star formation process that 
creates a stellar population with initial dynamical properties in¬ 
herited from their parent gas elements (SPH particles). During 
the first Gyrs, when the secondary bar grows, gas inflow and 
star formation are responsible for bringing energy to the inner 
waves. When the central star formation fades out as the mass of 
gas inflowing in the central region decreases, the energy dissi¬ 
pation at ILRs prevails. This is basically the scenario of SimS 
for t k 2 Gyr. What delays the inner bar extinction in SimL is a 
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Fig. 14. Top and middle panels: as for Fig. 1^ at t = 7.6 Gyr, with 
(SimL, top) and without (SimL""^^, middle) star formation. The surface 
brightness ranges from 22.8 to 17 Bmagpc^^ for SimL and 22.8 to 19.1 
Bmagpc^^ for SimL’'°^^. Bottom: m = 2 power spectra for SirnL"®^^, 
centred at f = 7.6 (top panel). This figure can be directly compared with 
Fig-O 


sustained local star formation rate for several Gyr (cf. Fig 
the same way as for sustaining the cr-drop phenomenon 


igm 

in 


In 


niak & Champavert|2006 1, the regular feeding of the central re¬ 
gion with recent stellar populations enables secondary bars to 
be long-lived. Simulations by |Friedli & Martinet| ( |1993] ) did not 
take star formation into account, leading to the extinction of the 
secondary bar in roughly 5 turns. This is also the case of other 
past simulations with a gaseous component that have been dis¬ 
cussed before. A noticeable exception is Friedli et al. ( 1996| l 
who ran a simulation with star formation, but only on a short 
timescale (less than 2 Gyr). 

To firmly establish the role of star formation, we recomputed 
SimL from t = 5.25 Gyr until f = 7.6 Gyr switching off the 
star formation process in the simulation code (SimL"°®^ in short). 
Wozniak & Champavert ( |2006| l used this method for cr-drops. 
Figure 14 shows the gas and stellar mass distributions at f = 
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7.6 Gyr. They are significantly different from SimL. The gas 
distribution is more concentrated in SimL™®^ than in SimL since 
the gas is no longer consumed by the formation of new stars. 
This has two immediate consequences. First, this predictably 
makes dissolving the inner bar easier as the gas can now reach 
the centre, increasing there the mass that in turn axisymmetrizes 
the gravitational potential. Second, no more gas is turned into 
stars so that the mass of the inner stellar bar cannot increase. 
After roughly 2 Gyr, the inner bar is indeed barely detectable in 
the stellar distribution of SimL"°®^. The spectrogram centred at 
t - 1.6 Gyr, as in Fig.[^ confirms that, although it has not yet 
fully disappeared, the amplitude of the m - 2 mode has strongly 
decreased with respect to SimL. 

The role of star formation in sustaining the nuclear bar over 
several Gyr is thus demonstrated for the first time. 


6. Conclusions 

Unlike the outcome of most numerical simulations with both 
stellar and gaseous components, we have successfully simulated 
a long-lived inner bar embedded in a large-scale primary bar. 

The ratio of the two bar lengths, the ratio of pattern speeds, as 
well as the age of the inner stellar bar population and the central 
gas mass fit well with observations published in the literature. 
Moreover, throughout the simulation, our models go through 
various morphological phases representative of the diversity of 
observations, including SBO. 

The most important difference with past simulations lead¬ 
ing to short-lived double-barred galaxies is the lack of overlap 
between the primary bar inner Lindblad resonance (ILRp) and 
the nuclear bar corotation (CRs) and, more generally, the lack of 
any kind of resonance overlap. The absence of mode coupling, 
confirmed by a Fourier analysis, implies that to sustain a perma¬ 
nent nuclear bar mode another physical or dynamical mechanism 
must feed the central waves. 

Star formation in the central region is identified as possibly 
being responsible for bringing energy to the nuclear mode. Star 
formation is also responsible for regulating the gas mass accu¬ 
mulation close to the centre, in the sense that it prevents a strong 
increase in mass density that can destabilize the inner bar orbits. 

As a direct consequence (whereas the inner bar can be tem¬ 
porarily undetectable, leading the pattern speed to be barely 
measurable or unmeasurable), the corresponding perturbation 
modes survive and allow the revival of the nuclear structure less 
than 1 Gyr after it disappears, provided that star formation con¬ 
tinues. 

A side result of our study is that an overlap between the ILR 
of the main bar and the CR of the nuclear bar cannot be sys¬ 
tematically assumed as it is not a necessary condition for the ex¬ 
istence of double-barred galaxies. Another direct consequence 
of our simulations is the evidence that several morphological 
features in barred galaxies that have been identified by various 
names (double-bars, pseudo-bulges, triaxial bulges, etc.) origi¬ 
nate dynamically from an unique inner mode. SimL displays all 
these features at various stages of its evolution. This speaks for 
a global analysis of triaxiality in the central regions of barred 
galaxies as being a unique dynamical phenomenon due to per¬ 
sistent modes in the central region. 
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